Finite-Temperature Behavior of an Inter-species Fermionic Superfluid with Population Imbalance 
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We determine the superfluid transition temperature Tc and related finite temperature phase diagrams for the 
entire BCS-Bose Einstein condensation crossover in a homogeneous mixture of ®Li and *'^K atoms with pop- 
ulation imbalance. Our work is motivated by the recent observation of an inter-species Feshbach resonance. 
Pairing fluctuation effects, which significantly reduce Tc from the onset temperature for pairing (T*), provide 
reasonable estimates of Tc and indicate that the inter-species superfluid phase should be accessible in future 
experiments. Although a generalized-Sarma phase is not stable in the ground state near unitarity, our phase 
diagrams show that it appears as an intermediate-temperature superfluid. 
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Ultracold Fermi gases with tunable attractive interactions 
provide an exciting opportunity to study superfluidity in the 
very general context of a crossover from BCS theory to Bose- 
Einstein condensation (EEC). While initial experiments ad- 
dressed the more conventional situation of intra-species pair- 
ing with equal populations of the two effective spin states, 
more recently there has been an emphasis on population im- 
balanced gases HI |2> S B S. Many interesting phases 
have been contemplated including two which appear to have 
been observed in the laboratory: the so-called Sarma phase 

01 which represents a homogeneously polarized state, and 
a phase separated (or inhomogeneously polarized) state L8]. 

Adding to the excitement is the possibility of discovering 
a new form of superfluid involving inter-species pairing. A 
first step en route to the discovery is the recent observation of 
Feshbach resonances between ^Li and ''"K atoms [9]. If the 
transition temperatures are accessible, this tunable attractive 
interaction should enable BCS-BEC crossover in superfluid 
phases associated with unequal mass pairing. Moreover, it is 
of particular interest for its relevance to color superconductiv- 
ity in quark matter HQ] . 

In this paper we determine the transition temperatures for 
the entire BCS-BEC crossover in a homogeneous mixture of 
^Li and ^''K atoms with population imbalance. In addition we 
address the temperature-polarization phase diagrams associ- 
ated with inter-species superfluid phases at and around unitar- 
ity. We consider only a homogeneously polarized superfluid 
and exclude from consideration the phase separated state prin- 
cipally because it is now clear jstl that the normal regions in 
this heterogeneous phase correspond to a complicated corre- 
lated normal state, which is currently difficult to characterize 
at an analytical level, except in the Umit of extreme popula- 
tion imbalance. This normal state, which is distinct from a 
free Fermi gas, has been addressed numerically using quan- 
tum Monte Carlo simulations jl ill . 

Our T ^ Q calculations are performed in a fashion consis- 
tent with a generalized BCS-Leggett ground state which has 



been studied previously 11121 11311 in the strict T = limit. 
These calculations show the importance of including pairing 
fluctuations which greatly suppress the transition temperature 
Tc from the pair formation temperature T*. An understand- 
ing of finite temperature effects positions us to address actual 
experiments (which are never in the ground state). Moreover, 
new superfluid phases appear which are constrained to an in- 
termediate regime of non-zero temperatures. That is, they are 
associated with a lower as well as upper critical temperature. 
We show that, in the absence of a trap, this intermediate tem- 
perature superfluid will be extremely difficult to observe when 
the heavy species is the majority, but it should be more acces- 
sible for the case where the heavy species is the minority. Fi- 
nally, we study how the phase diagram evolves as one crosses 
from BCS to BEC. In contrast to polarized Fermi gases with 
equal mass, close to, but on the BEC side of resonance, the in- 
termediate temperature superfluid disappears when the lighter 
species is the majority, giving way to a conventional polarized 
superfluid with only one transition temperature. 

There have been extensive studies on zero temperature 
properties of homogeneous fl^l as well as trapped inter- 
species Fermi gases [12, 14, 15]. This body of work (like that 
in the present paper) is based on a natural generalization of 
the BCS-Leggett ground state to accommodate unequal popu- 
lations. Similarly, while a two channel model formalism may 
be more relevant to the narrow resonances seen in Ref. |@1, 
all work to date (including our own) deals with the simpler 
one channel model. Theoretical studies at finite temperatures 
that are consistent with these T = calculations has been 
limited to a strict mean field approach which ignores the 
important effects of pair fluctuations, or non-condensed pairs. 

We begin with an outline of the central equations associ- 
ated with a generalized "Sarma" state appropriate to the case 
of unequal masses and finite polarization. This is followed 
by a short microscopic derivation. We choose the conven- 
tion > m|, so that the mass of spin-down fermions is the 
larger The mass ratio is mi/m-^ = 6.7. We define i?fe|,x = 



2 



Ek ± 4", Ek = + where = (6t ± a-i)/2- 

Here efc^- - /i^, and 6^^^ = P/2ma, where cr =t, J- 

The four unknowns which must be determined at general tem- 
perature T involve the two fermionic chemical potentials fi-^, 
Hi, and the excitation gap A, whose square appears in Ek- 
The quantity contains a contribution from condensed (sc) 
and non-condensed (pg) pairs: 



A^ 



(1) 



so that one of the gap components (Apg) must be separately 
determined in order to establish the transition temperature Tc- 
This is the lowest temperature at which A^c vanishes. 

There are then four equations, three of which can be associ- 
ated with strict mean field theory, and all of which are derived 
microscopically below. The equations for the total number 



n = 
are 



and 



and number difference Sn 



of fermions 



Ek- 



6n^J2^f{Eki)-f{Ek^)]. 



(3) 



Here f{x) — (e^^^ + 1) ^ is the Fermi distribution function. 
The gap parameter A is obtained from 



M 
2'Ka 



n- 



f{Eki)-f{Eki) 



2Eh 
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(4) 



Here the coupling constant is regularized by g 



where a is the s-wave scattering 
mi) is the reduced mass and 



Af/(27ra) - Ek(2ei 
length, M = m|m|/(m| 
ek = k^/2M. 

Finally, we need an equation for A^^ which requires that 
we establish the dispersion of the non-condensed pairs. These 
non-condensed pairs, or pseudogap effects, appear at T 7^ 
and are included via a T-matrix contribution to the fermion 
self-energy. Following Refs. It7k il61 . the fermionic self energy 
Ti„{K) = t{Q)Ga{Q — K), where the four-momentum 
Q = (if7;,q), K = (iw„,k), and Vli ( Un ) is the bosonic 
(fermionic) Matsubara frequency, with = T 

= ^X^nSk' ^ = The T-matrix is pre- 

sumed to have the structure t{Q) ~ tsc{Q) + tpg{Q). The 
condensate contribution satisfies tsc{Q) = —{Al^/T)S{Q). 
Here the fermionic Green's function is (K) = [Gq^ (K) — 
S,(i^)]-i, with Go^{K) = (zc^„ - 6.). We set h = 1 and 
kB = 1. 

The excited pair propagator (which is what we need to ar- 
rive at our fourth equation) is given hy tpg{Q) = [g^^ + 
x{Q)]~^, where the symmetrized pair susceptibility, xiQ) — 
J2k (7 ^OcriQ — K)Ga{K)/2, is used. A central assumption 
is the usual BEC condition that the pair chemical potential 
vanishes below Tc, which will, in turn, lead to Eq. This 
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Figure 1: (Color online) T* as a function of l/{kFa) for several 
values of p (as labelled). Here p < when the lighter (spin up) 
species is the majority (left panel) and p > otherwise. Insets: T* 
(Black solid line) and unstable regime (shaded regime) forp = —0.5 
andp = 0.5, respectively. 



BEC condition implies that tpg is dominated by terms with 
0. Importantly, the pseudogap is determined by Apg = 



J^Q tpgiQ)- It follows that T,„{K) = ~A^Gosi-K), and 
in that way we have derived Eq. ([T]i. One arrives at the two 
equations for the number densities via ~ J2k (^) ^^'^ 
in this way derives Eqs. (|2|i and (|3]l. 



The T-matrix may be expanded jlTIl as t~g{Q) — aoQ + 
aifl^ — £,^q^, after analytic continuation (ifli ^ + iO'^), 
where we have neglected the small imaginary part Fg. The 
pseudogap contribution can be written as 



E 



(5) 



Here b{x) is the Bose distribution function and (iq = 
(■\/aQ + 4ai^^(7^— ao)/2ai. In the BEC limit, it can be shown 
that ai/ao and flq q^/2M*, where M* = 
is the effective pair mass. Importantly, in this limit M* ap- 
proaches the total mass of the two constituent fermions, and 
Tc approaches the BEC temperature of ideal bosons of density 
min(n|, ni)/2 and mass M* . 

The gap equation, Eq. (|4|i, is equivalent to an extremal 
condition on the thermodynamic potential OVIm f / dA = 0, 
where 



n 



MF 



k,(T 



In 1 



(6) 

Superfluid stability requires that the number susceptibility 
matrix dn^r/dfia' should be positive-definite ff]. This can be 
shown to coincide with the condition that d^Vli^ip / dA^ > 0. 
When this condition is violated, phase separation may oc- 
cur These alternative phases have been widely studied ifisll 
at T = 0, as well as at non-zero temperature [ 19[]. 

The distinction between A and the order parameter Asc is 
an important component of the present theory. The former is 
associated with an onset temperature T* and the latter with 
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Figure 2: (Color online) Tc as a function of l/{kFa) for selected 
values of p. Here p < when the lighter species is the majority 
(left panel) and p > otherwise. The polarized superfluid solution 
is unstable in the shaded regions. 



Tc- We aiTive at a reasonable estimate of T* from Eqs. (|2]l- 
(HJi which is plotted in Fig. [T] as a function of 1/ (kpa). Two 
different signs of the polarization p — Sn/n aie indicated in 
the right and left panels. Here kp and Tp are the Fermi mo- 
mentum and the Fermi temperature of an unpolarized non- 
interacting Fermi gas with the same total particle density pre- 
suming a mass equal to the average mass of ^Li and ""^K. 
We will see that throughout the paper a superfluid phase with 
p < (the lighter species is the majority) appears to be more 
readily obtainable than one with p > 0. The arguments be- 
hind this asymmetry in p are subtle and involve both energetic 
comparisons as well as mechanical stability (in the sense of 
O^Umf /9A^ > 0). It can be analytically shown that at 
T = when polarization is small in the BEC regime, a su- 
perfluid phase with p < is more energetically favorable than 
one with p > 0. 

The figure shows that T* vanishes when the attraction is 
sufficiently weak and near its vanishing point it displays non- 
monotonic behavior Similar behavior has been observed 
previously for population imbalanced Fermi gases of equal 
masses |7]. The insets of Fig. [T] indicate the unstable regimes 
for p — ±0.5. The unstable regimes are asymmetric in the 
sign of p and at T = the behavior is consistent with phase 
diagrams obtained earlier lll2lll3ll . 

We turn next to the superfluid transition temperature, Tc, 
which is plotted in Fig.|2]as a function of l/(fci?a), for both 
p < (left panel) and p > (right panel). It is clear that 
these homogeneous superfluid phases are more likely to be 
observed when polarization is low. The shaded regions indi- 
cate where this form of superfluidity is unstable. Note that 
there is a rather pronounced asymmetry between the p < 
and p > cases. Indeed, when p > 0, a stable superfluid 
cannot be found near 1/kpa — 1.5 although it will emerge 
again deep in BEC. An "intermediate temperature superfluid" 
phase exists in the BCS through unitary regimes, which is sta- 
ble only away from the ground state. This is an unusual phase 
which has both a lower and an upper critical temperature, as 
previously found for population imbalanced same mass Fermi 




Figure 3: (Color online) Phase diagrams for mixtures of "Li and ''"K 
atoms at unitarity. Here p < when ^Li is the majority species 
and p > otherwise. The black solid, red solid, and green dashed 
lines are T* , Tc, and the boundary of stable phases, respectively. 
Labeled are polarized superfluid (brown), pseudogap state (yellow) 
and unstable Sarma (grey) phases. The white open space is polarized 
normal Fermi gas. 



gases 101 ■ 

Since the strongly interacting, unitary regime has been the 
focus of much interest, in Fig. [3] we first present the phase 
diagram at unitarity as a function of temperature and polar- 
ization p. It shows that uniformly polarized superfluid exists 
at low IpI and low but finite T (the dark shaded region). At 
higher T and higher |p|, there is non-superfluid state with a 
finite excitation gap (called the pseudogap) before the system 
becomes essentially uncoiTelated which we refer to as the nor- 
mal Fermi gas state. However, at low T and relatively high 
|p|, the Sarma superfluid phase (light shaded) is found to be 
unstable. Interestingly, when the heavier species is the major- 
ity (i.e., p > 0), the stable or unstable paired states occupy 
a smaller phase space than the opposite case. It should be 
noted that when phase separation is included, it may cover the 
entire unstable region and extend slightly into the pseudogap 
and superfluid regions in the phase diagram. In this way the 
re-entrant behavior (seen upon a vertical cut a constant low p) 
shown in the right hand panel of Fig. [3] should not survive. 

The boundary separating the pseudogap phase and the nor- 
mal Fermi gas can be associated with T*, while the line sepa- 
rating the polarized superfluid and the pseudogap phase corre- 
sponds to Tf.. In a strict mean field calculation [13] there is no 
pseudogap so that the T^. curve falls right on top of the T* line. 
It is clear, then, that pair fluctuation effects are extremely im- 
portant for they greatly reduce the regime where superfluidity 
appears. In the equal mass case, where there is an opportu- 
nity to compare with experiments lEoll . a similar reduction in 
the regime of stable superfluid was found iiHl which was in 
quite good agreement with the data. One can see from Fig. [3] 
that, for both p < and p > 0, while mean-field theory again 
predicts a large value of the upper critical polarization (be- 
yond which the superfluid phase cannot exist) and does not 
distinguish T* and Tc, our theory which includes pairing fluc- 
tuations predicts a smaller stable superfluid regime. Based 
on evidence in the equal-mass case, our prediction should be 
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over, we have previously ll2lll characterized the changes in the 
phase diagram for the equal mass case upon going from the 
homogeneous to the trapped situation. We have shown that, 
while the regimes of stability of generalized Sarma phases are 
greatly expanded in a trap, the characteristic values of the tran- 
sition temperatures are not significantly altered. 

This work is supported by NSF PHY-0555325 and NSF- 
MRSEC Grant 0820054. We thank Cheng Chin for helpful 
conversations. 



Figure 4: (Color online) Phase diagrams for mixtures of "Li and "'"K 
atoms at l/kpa = 0.5. The convention follows that in Fig.[3] 



more typical of future experiments. 

Another notable feature of the figure is that when the heav- 
ier species is the majority, a stable homogeneously polarized 
superfluid only exists in a very narrow regime of extremely 
low polarization. This should serve as an important guide to 
future experiments, for it suggests one has to be very careful 
in order not to miss the (homogeneously) polarized superfluid 
in this case. It might seem as though the transition tempera- 
tures are much higher in the case of unequal masses than in 
the equal mass case. We stress that it is more meaningful to 
compare quantities such as Tc/T* than Tc/Tp because, un- 
like in the equal mass case, here the energy unit Tp does not 
correspond to the Fermi energy of ether species of the atoms. 
For the equal-mass case Tc/T* « 0.5 while T^/T* « 0.3 for 
a mixture of ^Li and ''"K atoms at p = 0. 

To address how the phase diagram evolves from unitarity to 
the BEC side of the Feshbach resonance, we present in Fig.|4] 
the counterpart phase diagrams at l/kpa = 0.5 for the inter- 
species superfluid. It should be clear from the figure that, just 
as in the previous case, pair fluctuation effects which allow us 
(via Eq. ([T]i) to distinguish between the gap A and the order 
parameter Age, are extremely important as they greatly reduce 
the regime of stable homogeneous superfluidity. 

When the heavier species is the majority (p > 0), there is 
virtually no stable polarized superfluid. Again, the re-entrant 
behavior of the intermediate temperature superfluid at low 
p > in Fig. m should not survive if phase separation is con- 
sidered. In contrast, when the lighter species is the majority 
(p < 0), there is no longer an intermediate temperature su- 
perfluid phase. Rather a stable superfluid can be found for all 
polarizations and temperatures below the (single) critical tem- 
perature Tc- This behavior can be contrasted with the equal 
mass case where, for, e.g., p = 0.5 at T = 0, a stable su- 
perfluid cannot be found until deep in the BEC regime, when 
l/kpa > 2. For a moderate polarization, say p — —0.5, this 
transition temperature Tc « O.IST^, which should be acces- 
sible in future. 

Although we have considered the homogeneous rather than 
a trapped case, there are now experimental capabilities for ad- 
dressing this phase diagram using tomography ll22ll . More- 
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